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ABSTRACT 

We present the GalMer database, a library of galaxy merger simulations, made available to users through tools compatible with the 
Virtual Observatory (VO) standards adapted specially for this theoretical database. To investigate the physics of galaxy formation 
through hierarchical merging, it is necessary to simulate galaxy interactions varying a large number of parameters: morphological 
types, mass ratios, orbital configurations, etc. On one side, these simulations have to be run in a cosmological context, able to provide 
a large number of galaxy pairs, with boundary conditions given by the large-scale simulations, on the other side the resolution has 
to be high enough at galaxy scales, to provide realistic physics. The GalMer database is a library of thousands simulations of galaxy 
mergers at moderate spatial resolution and it is a compromise between the diversity of initial conditions and the details of underlying 
physics. We provide all coordinates and data of simulated particles in FITS binary tables. The main advantages of the database are VO 
access interfaces and value-added services which allow users to compare the results of the simulations directly to observations: stellar 
population modelling, dust extinction, spectra, images, visualisation using dedicated VO tools. The GalMer value-added services 
can be used as virtual telescope producing broadband images, ID spectra, 3D spectral datacubes, thus making our database oriented 
towards the usage by observers. We present several examples of the GalMer database scientific usage obtained from the analysis of 
simulations and modelling their stellar population properties, including: (1) studies of the star formation efficiency in interactions; 
(2) creation of old counter-rotating components; (3) reshaping metallicity profiles in elliptical galaxies; (4) orbital to internal angular 
momentum transfer; (5) reproducing observed colour bimodality of galaxies. 
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1. Introduction 

In the framework of the present cosmological paradigm, merg- 
ers and interactions are among the most important mechanisms 
govern ing galaxy formation and evolution. Spitzer & Baad3 
(Il951h proposed that collisions of late-type disc galaxies should 
produce early-type ones. A natural consequence of this phe- 
nomenon is th e morphology- density relation discovered three 
decades later ( |Dressleil[T980l) : dense regions of the Universe, 
where galaxy collisions are supposed to be more frequent, con- 
tain larger fractions of early-type galaxies than sparsely popu- 

lat ed areas. 

iToomre & Toomr3 (Il972h were the first to suggest that in- 
teractions 'tend to bring deep into a galaxy a fairly sudden 
supply of fresh fuel in the form of interstellar material. . .' 
The gas in the bar formed during the interaction loses its 
angular momentum due to the torques dCombes et al.l [T990l: 
iBames & Hernqui"s3ll996h . falls onto the galaxy centre possi- 
bly inducing strong bursts of star formation ( iMihos & Hernquisi 
1 19944 19961) and creating young compact central stellar compo- 



nent OS & Hemauistill994bh often observed in present-da y 
early-type galaxies (ISiFchenkol 2006t iKuntschner et al.l l2006h . 
Intense star formation episodes accompanied by supernova ex- 
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plosions, enrich the interstellar medium (ISM) with heavy ele- 
ments, consumed into the stars hence increasing the observed 
metal abundances of the stellar population. 

Large scale cosmological A^-bo dy simulations (e.g. 
ISpringelet aDl2006t lOcvirk et al1l2008h often lack the spatial 
resolution to trace in detail star formation and morphological 
transformation at galaxy scale. Therefore, usually they are 
complemented by semi-analytical prescriptions qualitatively 
accounting for phenomena strongly affecting galaxy evolution, 
such as star formation (e.g. lBlaizot et aIj|2004tlSomervilIe et al.l 
200^. However, the parameters of the semi-analytical models 
have to be chosen based on more detailed simulations of 
galaxy interactions. H igh resolution galaxy simulations (e.g. 
iBournaud et al.l l2b08h cannot be performed in large statistical 
numbers. A compromise has then to be done between statistics 
and resolution. This becomes one of the main motivations for 
studying large numbers of galaxy interactions by means of 
dedicated intermediate-resolution numerical simulations. 

Merger-induced star formation, as well as morphological 
transformation, strongly depends on the mass ratio of the 
interacting galaxies. Generally, the intensity of sta rburst de- 
creas es as the merger mass ratio increases (e.g. ICox et all 
l2008h . Equal mass mergers of disc galaxies (mass ratios be- 
low 4:1) usua lly result in early-typ e elliptical-Iike remnants 
(lToomrelll977l:lNaab & BurkertllBol while minor mergers (ra- 
tios above 10:1) do not destroy the progenitor's disc preserv- 
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ing its exponential mass distribution although making it thicker 
and dynamically hotter ('Ouinn et al." 1993': ' Walker et al.lll996t 
[Velazquez & White 1999; Bournaud et al. 2005]). A sequence of 
repeated minor mergers can form elliptical galaxies, with global 
morphological and kine matical properties sim ilar to that ob- 
served in real ellipticals jBournaud et alJlToOTh . 

Orbital parameters of the interaction and orientation of 
galaxies also strongly affect the process of merger, e.g. star for- 
mation efficiency on retrograde orbits is ge nerally higher than 
for direct encounters dPi Matteo et al.ll200'7l) . 

One needs to explore a large multi-dimensional parameter 
space (different initial morphologies related to gas content, or- 
bital configurations, mass ratios, etc.) by running thousands of 
simulations in order to fully understand the astrophysical conse- 
quences of galaxy interactions for the modern picture of galaxy 
evolution. 

The GalMer project, developed in the framework of the 
French national horizonQ collaboration, has the ambitious 
goal of providing access for the astronomical community to the 
results of massive intermediate-resolution numerical simulations 
of galaxy interactions in pairs, covering as much as possible the 
parameter space of the initial conditions and, thus allowing to 
study statistically the star formation enhancements, structural 
and dynamical properties of merger remnants. An important as- 
pect is to integrate the data services into the framework of the 
International Virtual Observatory in order to take the full ad- 
vantage of already developed technologies and data visualisation 
and processing tools. 

This paper is a technical presentation of: (1) the GalMer 
TreeSPH simulations providing all essential details on initial 
conditions for galaxies of diff'erent masses and prescriptions 
used to model the processes of star formation including su- 
pernova feedback and metallicity evolution; (2) the GalMer 
database, the first VO resource containing results of TreeSPH 
simulations; (3) the GalMer value-added services aimed at fa- 
cilitating the comparison between simulations and observations 
such as modelling the spectrophotometric galaxy properties us- 
ing an evolutionary synthesis code. 

The paper is organised as follows: in Section 2 we describe 
initial conditions of the numerical simulations and orbital pa- 
rameters of galaxy interactions; Section 3 contains information 
on the numerical method, prescriptions for star formation and 
metallicity evolution; Section 4 provides the description of the 
GalMer database structure, its access interface, and mechanisms 
of data visualisation; Section 5 presents value-added services of 
the GalMer database; in Section 6 we define some astrophys- 
ical applications which can be tackled with our simulations; 
Section 7 contains a brief summary. 

2. Initial conditions 

2.1. Galaxy models: moving along the Hubble sequence 

We model interactions among galaxies of different morpholo- 
gies, from ellipticals to late-type spirals. 

The adopted galaxy models consist of a spherical non- 
rotating dark-matter halo, which may or may not contain a stel- 
lar and a gaseous disc and, optionally, a central non-rotating 
bulge. For each galaxy type, the halo and the optional bulge 
are modelled as Plummer spheres, with characteristic masses 
Mb and Mh and characteristic radii and r//. Our choice of 
adopting a core density distribution for the dark halo seems 

' http://www.projet-horizon.fr/ 



to be more in accordance with the rotation curves of local 
spirals and dwarf galaxies, than the cuspy profiles predicted 
by Cold Dark Matter simulations (see |Di Matteo et al. 20083, 
Sec t. 2.4. 3 for a discussion). Th e stellar and gaseous discs follow 
the iMivamoto & Nagail (1 19751) density profile, with masses M, 
and Mg and vertical and radial scale lengths given, respectively, 
by /i, and a,, and hg and ag. We refer the reader to Appendix 
|A]for the analytical expression of these profiles. For the initial 
models of the disc galaxies, we chose an initial Toomre param- 
eter for the stellar disc Qstai = 1 -2, and two different initial ones 
for the gas component, Qgas = 0.3 and Qgas = 1.2, in order 
to reproduce, at least partially, the variety of parameters found 
in real galaxies jMartin & Kennicutil200ll : lBoissier et~al ] l200l 
iHitschfeld et al.ll2009h . ^ 

The database contains interacting galaxy pairs of different 
mass ratios (1:1, 1:2, 1:10), involving a giant galaxy (gEO for 
a giant-like elliptical, gSO for a giant-like lenticular, gSa for a 
giant-like Sa spiral, gSb for a giant-like Sbc spiral and gSd for a 
giant-like Sd spiral), interacting with: 

- either another giant galaxy; 

- or an intermediate-mass galaxy (hereafter iEO, iSO, iSa, iSb 
and iSd), having a total mass half that of the giant's mass; 

- or a dwarf galaxy (hereafter dEO, dSO, dSa, dSb and dSd) 
whose total mass is ten times smaller that that of the giant 
galaxy. 

The complete Hubble sequence of the galaxy models is 
given in FiglT] Moving along the Hubble sequence, from gEO 
to gSd galaxies, the mass of the central spheroid varies from 
/Wb = 1-6 X 10"/Wq for a gEO to Mb -0 for a gSd, while the gas 
mass Mg, absent in the case of a gEO and gSO, increases from 
9.2 X IO'^/Wq in a gSa to 1.7 x lO'^/Wg for a gSd (see Table [T] 
for a complete list of all the parameters of giant-like galaxies, 
and Tables|2]and[3]for the corresponding set of parameters of in- 
termediate and dwarf systems). Our modelled Hubble sequence 
capture the main properties of local galaxies. We plan in the fu- 
ture to extend this library to higher redshift, considering larger 
gas fractions for the disc galaxies. 

The initial rotation curves of the modelled disc galaxies are 
shown in Fig|2l 

Each galaxy is made up of a total number of particles /Vjot 
distributed among gas, stars and dark matter Gas particles are 
actually "hybrid particles", characterized by two mass values: 
one corresponds to the gravitational mass and stays unchanged 
during the whole simulation, and the other is the gas content of 
the particule, decreasing or increasing according to the local star 
formation rate and mass loss (see Sect. 3.1, Eq. 1). When the gas 
fraction is below a certain threshold, the particle is transformed 
into a star paiticle, and the remaining gas mass is distributed over 
neighbouring gas particles. 

For giant-giant interactions, each galaxy is made up of a 
total number of particles of /Vtot= 120000, distributed among 
gas, stars and dark matter, depending on the morphology (Table 
Ell. For giant-intermediate and giant-dwarf interactions, we in- 
creased by a factor 4 the total number of particles in the gi- 
ant galaxy. This allowed us to improve the spatial resolution 
of the simulations and, in particular, to maintain a high enough 
numerical resolution for the smaller galaxy (/Vtot=240000 and 
A^TOT=48000 for the intermediate and dwarf galaxy, respec- 
tively; Tables |5] and |6]l. 

To initialize particle velocities, we adopted the method de- 
scribed in Hernquist (1993). 

In order to distinguish the role of interactions from secular 
evolution, for each galaxy model we provide in the database 



Chilingarian et al.: The GalMer database 



1 ■ 1 ■ 1 

■1 ■~~ 

1.1.4 


1 ■ 1 ■ 1 

aso - 


1 ■ 1 ■ i 

asa - 


1 ■ 1 ■ 1 
est - 


1 ■ 1 ■ 1 

flSd - 


1 . 1 . 4 


- - 

I.I.I 





- 




-' ■ 


-' iso'- 

.i 

«. 




1 ■ 1 ■ 1 

ISb - 




"i . 1 ' . i' 


■l l' . i' 


i . i >• 1 


'l " .' 1 ■ ' . (■ 


1 ■ 1 ■ 1 
dEO ' 

: # 

1.1.4 


1 I 1 I 1 
d50 - 

: ■* ; 

1 . 1 . 4 


1 I 1 I . 

dSft - 

ii'rdiL''' 




1 ■ 1 ■ 1 

dSb - 

- * - 

1.1.4 


1 ■ 1 ■ 1 

■m 

1.1.4 



Fig. 1. Hubble sequence of the GalMer galaxy models. From left to right, and from the top to the bottom, projection on the x-z plane 
of giant, intermediate and dwarf galaxies. Different colors correspond to the different galaxy components: dark matter (grey), bulge 
(red), stellar disc (orange), and gaseous disc (blue). Each frame is 50 kpc x 50 kpc in size. 
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Fig. 2. Rotation curves for the different galaxy models. The ini- 
tial rotation curve is shown (solid thick black curve), together 
with the contribution by the bulge component (thin solid red 
curve), the stellar disc (orange dashed curve), the gaseous disc 
(dot-dashed blue curve) and the dark matter halo (dotted grey 
curve). Lengths are in kpc and velocities in units of 100 km/s. 



also its evolution in isolation, for 3 Gyr. A brief description 
about how the modeled galaxies evolve in isolation is given in 
Appendi?jBl 



Table 1. Galaxy parameters for giant galaxies. The bulge and 
the halo are modelled as Plummer spheres, with characteristic 
masses Mb and Mh and characteristic radii rg and rn- M, and 
Mg represent the masses of the stellar and gaseous discs whose 
vertical and radial scale lengths are given, respectively, by /j, and 
a*, and hg and Og. 





gEO 


gso 


gSa 


gSb 


gSd 


Mb [2.3 X 10'' Mo] 


70 


10 


10 


5 





Mh [2.3 X lO^Mo] 


30 


50 


50 


75 


75 


M, [2.3 X 10' Mq] 





40 


40 


20 


25 


Mg/M. 






0.1 


0.2 


0.3 


ra [kpc] 


4 


2 


2 


1 




rn [kpc] 


7 


10 


10 


12 


15 


a, [kpc] 




4 


4 


5 


6 


h, [kpc] 




0.5 


0.5 


0.5 


0.5 


cig [kpc] 






5 


6 


7 


hg [kpc] 






0.2 


0.2 


0.2 



2.2. Orbital parameters 

For each pair of interacting galaxies, we performed several sim- 
ulations, varying the galaxies' orbital initial conditions (initial 
orbital energy E and angular momentum L) and taking into ac- 
count both direct and retrograde orbits (Tables |7][1] and |9]l. For 
each interacting pair, we kept the disc (when present) of one 
of the two galaxies in the orbital plane (/'i = 0°), and varied 
the inclination 12 of the companion disc, considering: 12 = 0°, 
i2 = 45°, 12 - 75°, and (2 - 90°. The clustering of the angles 
toward 12 - 90°, for an uneven sampling, is logical from a pure 



4 



Chilingarian et al.: The GalMer database 



Table 2. Galaxy parameters for intermediate galaxies. The bulge 
and the halo are modelled as Plummer spheres, with characteris- 
tic masses Mg and Mh and characteristic radii and r//. M« and 
Mg represent the masses of the stellar and gaseous discs whose 
vertical and radial scale lengths are given, respectively, by /j, and 
a,, and h„ and a„. 





iEO 


ISO 


iSa 


iSb 


iSd 


Mb [2.3 X 10" Mq] 


35 


5 


5 


2.5 





M« [2.3 X 10' Mq] 


15 


25 


25 


37.5 


37.5 


M. [2.3 X 10' Mo] 





20 


20 


10 


12.5 


MJM. 






0.1 


0.2 


0.3 


rg [kpc] 


2.8 


1.4 


1.4 


0.7 




rn [kpc] 


5. 


7. 


7. 


8.5 


10.6 


a, [kpc] 




2.8 


2.8 


3.5 


4.2 


h, [kpc] 




0.35 


0.35 


0.35 


0.35 


[kpc] 






3.5 


4.2 


5. 


hg [kpc] 






0.14 


0.14 


0.14 



Table 3. Galaxy parameters for dwarf galaxies. The bulge and 
the halo are modelled as Plummer spheres, with characteristic 
masses Mb and Mh and characteristic radii and rn- M,. and 
Mg represent the masses of the stellar and gaseous discs, whose 
vertical and radial scale lengths are given, respectively, by /j, and 



a*, and hg and ag. 





dEO 


dSO 


dSa 


dSb 


dSd 


Mb [2.3 X 10' Mo] 


7 


1 


1 


0.5 





Mh [2.3 X 10' Mo] 


3 


5 


5 


7.5 


7.5 


M. [2.3 X 10' Mo] 





4 


4 


2 


2.5 


Mg/M. 






0.1 


0.2 


0.3 


rg [kpc] 


1.3 


0.6 


0.6 


0.3 




r-H [kpc] 


2.2 


3.2 


3.2 


3.8 


4.7 


a, [kpc] 




1.3 


1.3 


1.6 


1.9 


ht [kpc] 




0.16 


0.16 


0.16 


0.16 


flg [kpc] 






1.6 


1.9 


2.2 


hg [kpc] 






0.06 


0.06 


0.06 



Table 4. Particle numbers for giant-giant interactions (mass ratio 
1:1) 





gEO 


gSO 


gSa 


gSb 


gSd 








20000 


40000 


60000 


A's.ar 


80000 


80000 


60000 


40000 


20000 




40000 


40000 


40000 


40000 


40000 



geometrical point of view, considering that the probability of the 
spin ;2 of the second galaxy to be oriented between and 12 is 
proportional to 1 - cos(/2): this means, for example, that an ori- 
entation i2 between 45° and 90° has probability 2.3 times higher 
than an orientation 12 between 0° and 45°. However, the spins 
aUgnment may not t) e totally uncorrelated, as recently shown by 
iJimenez et alj ( 1201 Oh . but a distribution function, to our knowl- 
edge, is still lacking. For giant-dwarf interactions, we also con- 
sider a more generic case, with ii - 33°, and i2 - 130° (see Fig. 
[3]for a sketch of the initial orbital geometry and Table [TOl for the 
orientation of the galaxy spins). 



3. Numerical method 

To model galaxy evolution, we employed a Tree-SPH code, in 
which gravitational forces are calculated using a hierarchical tree 



Table 5. Particle numbers for giant-intermediate interactions 

(mass ratio 1:2) 



gEO 



gSO 



gSa 



gSb 



gSd 



iVga.s 

Af„a,. 320000 
A^DM 160000 



320000 
160000 



80000 
240000 
160000 



160000 240000 
160000 80000 
160000 160000 



iEO 



iSO 



iSa 



iSb 



iSd 



A'star 160000 

A'dm 80000 



160000 
80000 



40000 
120000 
80000 



80000 120000 
80000 40000 
80000 80000 



Table 6. Particle numbers for giant-dwarf interactions (mass ra- 
tio 1:10) 



gEO 



gSO 



gSa 



gSb 



gSd 



A'sta, 320000 
A'dm 160000 



320000 
160000 



80000 
240000 
160000 



160000 240000 
160000 80000 
160000 160000 



dEO 



dSO 



dSa 



dSb 



dSd 



A'ga.s 

A^star 32000 
A'dm 16000 



8000 

32000 24000 
16000 16000 



16000 24000 
16000 8000 
16000 16000 




Fig. 3. Adopted orbital geometry for the simulations. We set up 
the collision in such a way that the orbital angular momentum 
is parallel to the z-axis and that the centers of the two galaxies 
initially are on the x-axis. The galaxy spins are represented by 
the blue and red arrows, respectively. They are specified in terms 
of the spherical coordinates (/i,<J>i) and (i2,<t2). See TablefTOlfor 
their initial values. 



method (iBarnes &Hutlll986h and gas evolution is followed by 
means of smoothed particle hydrodynamics (SPH, iLucvl [19771 : 
Gingold & Monaghan 1982). Gravitational forces are calculated 
using a tolerance parameter 9 - 0.7 and include terms up to the 
quadrupole order in the multipole expansion. A Plummer po- 
tential is used to soften gravitational forces, with same soften- 
ing lengths for all particle types. We assume a softening length 
e = 280 pc for giant-giant interactions, and e - 200 pc for the 
giant-intermediate and giant-dwarf runs. 

The code (evaluation of the gravitational forces, implemen- 
tation of the S PH technique and s t ar for mation modelling) 
was described in ISemelin & CombesI (l2002l) . where a standard 
validation test for this type of codes (the collapse of an ini- 
tially s tatic, isot hermal sphere of self-gra v itating gas, seelEvrardI 
(11988') and also'Hernauist & Katzl(ll989h : lThacker et al.1 (l2000l) . 
[Springel et al. 2001) was also presented. In our present study, we 
adopted an isot hermal gas phase for the G allVIer runs. Other tests 
are presented in lDi IVIatteo et al.l (l2008ah . 
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Table 7. Orbital parameters for giant-giant interactions Table 9. Orbital parameters for giant-dwarf interactions 



orb. id 






L= 




spin'^ 




kpc 


10-kms-' 


10^ km s^'kpc 


lO-^km^s-' 




Oldir 


100. 


2. 


56.6 


0.0 


up 


01 ret 


100. 


2. 


56.6 


0.0 


down 


02dir 


100. 


3. 


59.3 


2.5 


up 


02ret 


100. 


3. 


59.3 


2.5 


down 


03dir 


100. 


3.7 


62.0 


5.0 


up 


03ret 


100. 


3.7 


62.0 


5.0 


down 


04dir 


100. 


5.8 


71.5 


15.0 


up 


04ret 


100. 


5.8 


71.5 


15.0 


down 


05dir 


100. 


2. 


80.0 


0.0 


up 


05ret 


100. 


2. 


80.0 


0.0 


down 


06dir 


100. 


3. 


87.6 


2.5 


up 


06ret 


100. 


3. 


87.6 


2.5 


down 


07dir 


100. 


3.7 


94.6 


5.0 


up 


07ret 


100. 


3.7 


94.6 


5.0 


down 


OSdir 


100. 


5.8 


118.6 


15.0 


up 


08ret 


100. 


5.8 


118.6 


15.0 


down 


09dir 


100. 


2.0 


97.9 


0.0 


up 


09ret 


100. 


2.0 


97.9 


0.0 


down 


lOdir 


100. 


3.0 


111.7 


2.5 


up 


lOret 


100. 


3.0 


111.7 


2.5 


down 


lldir 


100. 


3.7 


123.9 


5.0 


up 


11 ret 


100. 


3.7 


123.9 


5.0 


down 


12dir 


100. 


5.8 


163.9 


15.0 


up 


12ret 


100. 


5.8 


163.9 


15.0 


down 



^ Initial distance between the two galaxies. 
Absolute value of the initial relative velocity. 

" L =1 Tini X Vim I- 

E = ViniV2 - G(mi -I- m2)/r,„i, with m, = mj = 2.3 x 10"Mo. 
° Orbital spin, if the z-component is parallel (up) or anti-parallel 
(down) to the z-axis 



orb.id 


''ini 






E'' 


spin"^ 




kpc 


lO^kms-' 


10^ km s^'kpc 


lO^km^s-^ 




Oldir 


100. 


1.48 


29.66 


0. 


up 


Olret 


100. 


1.48 


29.66 


0. 


down 


02dir 


100. 


1.52 


29.69 


0.05 


up 


02ret 


100. 


1.52 


29.69 


0.05 


down 


03dir 


100. 


1.55 


29.72 


0.1 


up 


03ret 


100. 


1.55 


29.72 


0.1 


down 


04dir 


100. 


1.48 


36.33 


0.0 


up 


04ret 


100. 


1.48 


36.33 


0.0 


down 


05dir 


100. 


1.52 


36.38 


0.05 


up 


05ret 


100. 


1.52 


36.38 


0.05 


down 


06dir 


100. 


1.55 


36.43 


0.1 


up 


06ret 


100. 


1.55 


36.43 


0.1 


down 



Initial distance between the two galaxies. 
*" Absolute value of the initial relative velocity. 
= L =1 ri„, X v,„i |. 

E = Vi„i2/2 - Gimi + m2)/nn., with mi = 2.3 x 10" Mq and mi = 

2.3 X IO'^Mq. 

Orbital spin, if the z-component is parallel (up) or anti-parallel 
(down) to the z-axis 



Table 10. Orientation of the galaxy spins, for giant-giant, giant- 
intermediate, and giant-dwarf simulations. 



giant-giant giant-dwarf 
giant-intermediate 
giant-dwarf 

T; W 33^ 

Oi 0° 0° 

(2 0°,45°,75°,90° 130° 

$2 0° 0° 



Smoothed Particle Hydrodynamics (SPH) is a Lagrangian 
technique in which the gas is partitioned into fluid elements rep- 
resented by particles, which obey equations of motion similar to 
the colUsionless component, but with additional terms describ- 



Table 8. Orbital parameters for giant-intermediate interactions 



orb.id 


ml 




L= 


£■1 


spin'^ 




kpc 


102 km g-i 


lO^km s^'kpc 






Oldir 


100. 


1.73 


48.0 


0.0 


up 


Olret 


100. 


1.73 


48.0 


0.0 


down 


02dir 


100. 


1.79 


49.0 


0.1 


up 


02ret 


100. 


1.79 


49.0 


0.1 


down 


03dir 


100. 


1.73 


60.0 


0.0 


up 


03ret 


100. 


1.73 


60.0 


0.0 


down 


04dir 


100. 


1.79 


61.0 


0.1 


up 


04ret 


100. 


1.79 


61.0 


0.1 


down 


05dir 


100. 


1.73 


69.3 


0.0 


up 


05ret 


100. 


1.73 


69.3 


0.0 


down 


06dir 


100. 


1.79 


69.7 


0.1 


up 


06ret 


100. 


1.79 


69.7 


0.1 


down 



Initial distance between the two galaxies. 
Absolute value of the initial relative velocity. 

L =1 Tmi X Vim I- 

E = Vim2/2 - G(mi + m2)/rim, with mi = 2.3 x 10" Mq and mj = 
1.15 x lO'^Mo. 

Orbital spin, if the z-component is parallel (up) or anti-parallel 
(down) to the z-axis 



ing pressure gradients, viscous forces and radiative effects in gas. 
To capture shocks a conventional form of the artificial viscosity 
is use d, with parameters ff = 0.5 and p - l.O dHernquist & Kata 
Ll989). To describe different spatial dynamical ranges, SPH par- 
ticles have individual smoothing lengths /i,, calculated in such 
a way that a constant number of neighbours Ues within 2/z,. 
The giant-giant simulations were performed using a number of 
neighbours A^^ ~ 15; for the giant-intermediate and giant-dwarf 
interactions A^, ~ 50. The gas is modelled as isothermal, with 
a temperature Tgas = IQi^K. Because of the short cooling time 
of disc gas, fluctuations in the gas temperature are quickly radi- 
ated away, so that simulations employing an isothermal equation 
of state differ little from more realistic ones (Mihos & HernquisJ 
119961: lNaabet"ani2006h . 

The equations of motion are integrated using a leapfrog al- 
gorithm with a fixed time step of Af = 5 x 10^ yr. 

3.1. Star Formation and continuous stellar mass loss 



Numerous prescriptions and 

19921: ISteinmetz & Muep ^ 



tech niques 
1991 



iw/t ISteinmetz & Muelle fl liv^4t lap: 
Springel & Hernquisli 120031: ICox etalT 1200^ 



exis t 
Springell 



for modelling 
star form ation rate and feedback in nu merical simulations. 

As in lMihos & Hernquisj(ll994bh . we parametrized the star 
formation efficiency for a SPH particle as 



-CXpgas 

pas 



(1) 
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with the constant C - 0.3 pc-'^^Mo^'^-Gyr"' such that the iso- 
lated giant disc galaxies form stars at an average rate of between 
1 and 2.5MQyr\ 

The choice of the parametrization in Eq.[T]is consistent with 
the observational evidence that on global scales the SFR in disc 
galaxies is well represented by a Schmidt law of the form Ssfr - 
ASgas'^, Sgas and Zsfr being disc-averaged surface densities, 
with the best fitting slope about 1.4 (see Kennicutt 1998,, but 
also IWong & Blitj 120021: iBoissier et al.l 120031: iGao & Solomon! 
I2004ir This relation see ms to apply, with a simi lar slope, also to 
local scales, as shown in lKennicutt et al.l (l2005h for Messier 51. 
As we checked, the prescriptions we adopted well reproduces 
the Schmidt law in t he Z^as - ^.sfr plane, both for isolated and 
interacting galaxies dPi Matteo et al.ll2007h . 

Once the SFR recipe is defined, we ap ply it to SPH parti- 
cles, us ing the hybrid method described in iMihos & HemquistI 
(Il994bh . In this method, each initial gas particle is in fact hy- 
brid containing a gas fraction and a star fraction. Its gravitational 
mass Mi remains constant, but the gas content of the particle, 
^/,gas, changes overtime according to Eq.[T] Gravitational forces 
are always evaluated using the gravitational mass M,, while hy- 
drodynamical quantities, in turn, use the time-varying M,;gas. If 
the gas fraction present in the hybrid particles drops below 5% 
of the initial gas content, the hybrid particle is totally converted 
into a star-like particle and the small amount of gas material still 
present is spread over the gas fraction of the neighbours. 

The effects of star formation on the surrounding 
ISM were implem e nted u sing the technique described in 
iMihos & Hernquisi (Il994bl) . For each star-forming hybrid 
particle, we evaluate the number of stars formed with masses 
> 8Mq (a [Miller & Scalo (1979ji) stellar initial mass function 
(IMF) is adopted), and we assume that they instantaneously 
become supernovae, leaving behind remnants of I.4M0 and 
releasing their mass to the surrounding ISM. The released mass 
also enriches the metallicity of the surrounding gas. This is 
done assuming a yield y = M,.f,/M, =0.02, where Mr^, is the 
total mass of all reprocessed metals and M« the total mass in 
stars. For each gas particle, mass and metals return is applied to 
the ;-th neighbour gas particle, using a weight w,- based on the 
smoothing kernel. 

The energy injection in the ISM from SNe explosions is ac- 
counted assuming that only a fraction ekin of Esn - 10^ 'erg 
goes into kinetic energy, by applying a radial kick to velocities 
of neighbour gas particles; thus, for each SNe explosion, the /- 
th neighbouring gas particle receives a velocity impulse directed 
radially away from the "donor" 



' 2 Wj eidn Esi< 
, Mi 



1/2 



(2) 



w, being the weighting based on the smoothing kernel and M, 
the mass of the receiver. 

The value of eyn is chosen such that the total amount of ki- 
netic energy received by a gas particle, due to the contribution 
from all its neighbours, is < 1 km s"' to prevent a rapid growth 
of the vertical thickness of the gaseous disc. 

Together with star formation, we also take into account the 
competing process of stellar mass-loss: at each time step, an 
amount 



MiAt) = 



(M,-M,-,gas(f))Af 



Co 



t - fbirth + Tq 



(3) 



the particle. In the formula above, fbinh represents the birth time 
of the population, Tq - 4.97 Myr and cq - 5.47 x 10"^ (see 
Jungwiert et al..,2001i, for details). 

3.2. Metallicity evolution 

The metallicity content of the modelled giant galaxies is initially 
distributed according to a radial profile 

z„,(R) = zo X l0-»'"*« (4) 

where R is the particle distance from the galaxy center and zo = 
3 xza (iKennicutt et al.ll2003HMagrini et al.ll2007HLemasle et all 
l2008h . For intermediate and dwarf galaxies, we generalised this 
formula, taken acco unt a dependency on the galaxy mass a nd on 
its half-mass radius (iTremonti et al.ll2004l: iLee et al.ll2006h : 



(5) 



of the stellar mass formed in the hybrid particle is assumed lost 
by evolutionary effects, going to enrich the gas content M,,gas of 



Z,„(R) = ^Mgal/Mgia„,Z0 X IQ-O-OT^-^^^R/ rso 
= ^Mgal/MgiantZo X lO-^-^^'^^''^" 

Mgai being the mass of the intermediate (or dwarf galaxy), rgo its 
half-mass radius, Mgiant the mass of the giant galaxy, and r^o = 
4.85 kpc the average half-mass radius of giant galaxies in the 
sample. 

The metallicity of the old stellar component is kept un- 
changed during the simulations, so that only remixing and dy- 
namical eff'ects can reshape the initial gradient of the old stellar 
population. In turn, the metallicity of the gas component and of 
the new stellar population (that formed during the simulation) 
changes with time, due to the release of metals from SNe explo- 
sions, as star formation proceeds. In more detail, hybrid particles 
are characterised by two metallicity values: z,„ and Znew The first 
value corresponds to the metallicity of the gas mass of the par- 
ticle, while the second one provides the metallicity of the stellar 
mass contained in the hybrid particle. As for old stars, initially 
their metallicity z,„ is distributed into the galaxy disc accordingly 
to Eq. m with the central regions more metal-rich than the outer 
disc. The metallicity Znew of a new stellar component, is set equal 
to that of the gas in which it forms, while the reprocessed met- 
als enrich the surrounding gas, according to the yield described 
before. 

4. The GalMer Database 

All the results of the GalMer simulations are accessible on- 
line using Virtual Observatory technologies. The three essential 
blocks for providing online access to the data are: (I) storing the 
data; (2) storing and querying the data description, i.e. metadata; 
(3) mechanisms for accessing and visualising the data. Besides, 
we provide services to perform online data analysis, which are 
described in detail in the next section. 



4. 1. Data format and storage 

Aimed at interopera bility and performance, we chose the FI TS 
Binary Table format (ICotton et alJll995HHanisch et alJl200lh to 
store the simulation results. The FITS format is handled by a 
variety of tools widely used by the astronomical community. A 
FITS bi nary table can be easily incorporated into the VOTable 
format dOchsenbein et al.ll2004t) used by new generation astro- 
nomical software tools introduced by the Virtual Observatory. 

Every galaxy interaction ("GalMer experiment" hereafter) 
includes 50 to 70 snapshots with a 50 Myr time interval con- 
taining data for individual particles traced by the simulations. 
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thus following the evolution of an interacting galaxy pair for 
2.5-3.5 Gyr We store each snapshot as an individual file. The 
following properties are provided for every particle: Cartesian 
coordinates {X, Y, Z), three-dimensional velocity vector (vx, vy, 
vz), total mass (M), particle type (hybrid, star or dark matter), 
identification of a galaxy which a given particle belonged to in 
the initial step of the simulation. Besides, we provide metallici- 
ties and fbiith, i-C- the average birth time of the stellar material in 
a given particle; for stellar particles, they are kept fixed through 
the simulations, since these particles do not evolve. For hybrid 
particles we provide mean fbinh and metallicities Zm instead, as 
well as gas masses, metallicities Znew of stars formed during the 
previous timestep (i.e. current gas metallicity). 

All the information related to the input parameters of a given 
GalMer experiment: morphological types of galaxies, orbital 
configuration, units of masses, coordinates, and velocities, as 
well as the epoch of a current snapshot are provided in FITS 
headers of snapshot files, making them self-consistent and avail- 
able for further stand-alone analysis without need to connect to 
the GalMer database. 

A typical snapshot containing 240000 particles has a size 
of 12 MBytes, resulting in 0.6-0.9 GBytes per GalMer exper- 
iment. The present data release contains simulations of about 
a thousand giant-giant interactions having a total volume of 
~0.9 TBytes. New simulations will be ingested into the database 
and put online as soon as they have been completed. 

4.2. Simulation metadata 

The metadata are computed for every individual snapshot at the 
time of the database update. Since the data are archived and do 
not change in time, neither the metadata do, the access to the 
database is read-only, unless new simulation results are ingested 
into it. 

The metadata of GalMer simulations conform to the current 
version of the SimDB data model (Lemson et al. in prep.) be- 
ing presently developed by the International Virtual Observatory 
Alliance (IVOA). SimDB is supposed to provide a complete self- 
sufficient description of A^-body simulations results using object- 
oriented approach and is designed using Unified Modeling 
Language (UML). For practical usage, e.g. for constructing a 
database containing the numerical simulation metadata, a UML 
data model has to be serialized. We partly serialize the SimDB 
data model into a relational database schema keeping another 
part (modified Characterisation class of SimDB) serialized into 
native XML. 

We use the open-source object-relational database man- 
agement system (DBMS) PostgreSQI^ to implement the ad- 
vanced metadata queryi ng mechanisms described in detail in 
IZolotukhin et al.l ( l2007h . We modified the original SimDB 
data model by replacing its Characterisation object with the 
full IVO A Characterisation Data Model (CharDM) metadata 
dLouvs e t al. 2008). CharDM is a way to say where, how ex- 
tended and in which way the observational or simulated dataset 
can be described in a multidimensional parameter space. Our 
metadata querying approach allows us to use additional WHERE 
clauses in standard SQL queries expressed as XPath statements 
to constrain particular CharDM elements without need to serial- 
ize rather complex CharDM structure into a relational database 
schema and, therefore, prevents additional complications of SQL 
queries. 



^ http://www.postgresql.org/ 



The relational DBMS stores metadata including links to 
FITS files containing all coordinates and data of individual sim- 
ulation particles. Therefore, to execute operations involving par- 
ticles, such as cutouts or statistics, it is necessary to fetch the ac- 
tual datasets. However, all global properties of simulations, such 
as total masses of gas and stars are available inside the DBMS 
making possible to extract a global star formation history of a 
given GalMer experiment at the database level. This function is 
implemented as a stored procedure inside the SQL database and 
is accessible from the web-interface. 

All actual metadata queries are executed by the server-side 
of the database web-interface using the parameters visualised in 
web-pages, therefore user does not need to type in SQL queries 
inside web-based forms. 

4.3. Data access and visualisation 

The interactive access to the data is provided through the WWW 
interface at http://galmerobspm.fr/. It is based on a de facto 
standard asynchronous JavaScript and XML {AJAX) technology 
and, thus, supports most widely-used modern Internet-browsers: 
Mozilla Firefox ver.>1.5, Microsoft Internet Explorer ver.>5, 
Apple Safari ver>3. Individual snapshots for every simulation 
can be directly accessed in a batch mode as well. 

The web-site provides a database query interface for access- 
ing simulations for given galaxy morphological types and or- 
bital parameters of interactions. Being a part of SimDB meta- 
data, this information is stored in the database, and is retrieved 
dynamically and displayed in a pop-up window, once user has 
selected particular elements in the query form. In order to use 
all available features of the GalMer database web-interface, the 
following software components have to be installed/enabled on 
the user's cornputer: (1) JavaScript support in a web browser; 
(2) SUN Jav(Q including the Java Applet browser plug-in and 
support for the Java WebStart funcionality which are normally 
configured automatically during the installation of SUN Java. 

The interactive data access includes several steps (see Fig.|4]i. 
At first, galaxy morphological types and orbital configurations 
of an interaction are chosen using the "DB Query" tab. Then, 
the user is asked to select one GalMer experiment in the "Query 
Results" tab from a list of those matching the selection criteria 
provided at the first step. Once an experiment has been selected, 
the user can download or visualise its integrated star formation 
history and/or download individual snapshots provided in the 
"Experiment" tab. Then, it is possible to preview the contents of 
a given snapshot and get access to the value-added data analysis 
tools for it using the "Shapshot" tab. 

In the "Snapshot" tab we provide a powerful AMZ-based 
preview mechanism for interactive data manipulation directly 
from within the web-browser (see bottom right panel in Fig.|4]i. 
The main purpose of the snapshot visualisation and manipulation 
is an interactive choice of the projection and scaling parameters 
for the generation of projected maps and synthetic images and/or 
spectra as described in the next section. 

Although the interactive data visualisation capabilities avail- 
able inside the web-browser remain limited, this limitation can 
be overpassed by using dedicated software tools if the data are 
sent to them directly from the database web-interface. 

We implemented the interaction between the GalMer 
database web-interface and existing Virtual Observatory tools 
dedicated for dealing with tables (topca^ lTavlodl2005h . im- 

^ http://www.java.com/ 

* http ://w w w. star. bri s . ac . uk/ ~mbt/topcat/ 
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Fig. 4. Interactive data access using the GalMer web-interface: (1) selecting galaxy morphologies and orbital types; (2) choosing 
one GalMer experiment from the list of those matching the selection criteria; (3) selecting a snapshot; (4) previewing the snapshot 
and accessing data analysis tools for it. 



ages (CDS AladiiB, 'Bonn arel etaP l2000h . and spectra (ESA 
VOSpe^ lOsuna et al.i ,2001" The mechanism of interaction 
was proposed bv IChilingarian & ZolotukhinI (l2008h and ini- 
tially implemented for the ob servational data archive ASPID- 
SR (Ch ilingarian et al. l20 7 [). It is based on the middleware 
dZolotuk hin & Chilingarianl BOOSl) represented by the two com- 
ponents (Fig. |5]l interacting using Sun Java LiveConnect: (1) 
a part handling user's actions inside the web-browser imple- 
mented in JavaScript and (2) a Java applet for sending data to 
the VO applications, implementing a VO application messag- 
ing protocol. Present version of the middleware implements the 
PLASTIC ( PLatform for Astronomical Tool Interconnection, 
iBoch et aP 12006) protocol. Support for the mo re advanced 
SAMP (Simple Application Messaging Protocol, iTavlor et"al] 
l2009h will be added in the next release of the database. 



5. Data Analysis Services 

We developed an application programming interface (API) li- 
brary to access and analyse the results of GalMer simulations. 
The mergerapi library is implemented in ANSI C and needs a 
minimal set of prerequisites to compile: the cfitsio FITS API 
librarjQ to access the simulation data and the GD librarjQ to 
generate gif and png output images directly viewable inside a 
web-browser The mergerapi library is installed on the server 
side and is used by several on-the-fly data analysis services de- 
scribed hereafter 

5.1. Projected Maps 

For a given snapshot we provide a set of services for the on-the- 
fly server-side computation of projected maps of various quan- 



^ http://aladin.u-strasbg.fr/ 
* http://esavo.esa.int/vospec/ 



' http://heasarc.gsfc.nasa.gov/docs/software/fitsio/fitsio.html 
* http://www.libgd.org/ 
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Middleware 
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<apiplet '3dc-''oirg.if[i-pra]ect.Ap3"orit'DLrtpp'^ct' nflin=-"flp3_ tout rol_ applet 
archiire^'/filej/app_caTtrol_app"let. jar'' 'rfidth^"3" heiclit^"!';' 



'a3p_co nt rcl_applei" 



JavaScript 

v/o_script.is 



Java 

-T,^ app_control_applet 



1. Controls app_controLapplet 1. isAladinStartedQ 

2. Launches VO client tools 2. isTopcatStartedO 

3. Checks permanently rf tools are avalable 3. IsVospecStartedQ 

4. Sends data to VO applications 4- sendSchptToAladinQ 

5. ... 



Virtual Observatory tools 





Fig. 5. The middleware connecting the GalMer database web- 
interface with the VO tools dedicated for advanced data manip- 
ulation at top, is illustrated a projection with mergerapi, and at 
bottom, a VOSPEC spectrum, and a map in Aladin. 



tities traced by the simulations. These tools are accessible at the 
"Snapshot" tab at the GalMer web-site. 

The following quantities can be computed for any combina- 
tion of the three particle types: surface density, line-of-sight ra- 
dial velocity (v,), line-of-sight velocity dispersion (cr). Gas mass, 
stellar mass, gas metallicity, and metallicity of new stars apply 
to hybrid particles only. Stellar metallicity provides a metallic- 
ity of underlying ("old") stellar population, therefore it applies 
to "star" particles and remains constant in time for a given par- 
ticle. Preview surface density maps are shown directly in a web- 
browser, while for other quantities (velocities, metallicities, etc.) 
the maps are available for download as FITS images or can be 
visualised with CDS Aladin. 

In order to generate maps on-the-fly, we developed an effi- 
cient computational algorithm. All the maps are computed for 
a parallel projection onto a plane (i.e. assuming an observer at 
the infinite distance) at the same time allowing to specify a spa- 
tial scale (i.e. to give an effective distance from an observer to 
the barycentre r). Therefore, the viewport is uniquely defined by 
the two quantities: azimuthal (<f) and polar {ff) angles. To sim- 
plify the comparison of computed maps with observations and 
mosaicing of maps with different spatial sampling (e.g. zoom-in 
on a given regions of a merger rem nant with the overall image) 
we use the FITS WCS convention dGreisen & Calabrettall2002h 
assuming a tangential projection and assigning right ascension 



and declination of zero to the projected barycentre position. In 
this case, the coordinates of a particle on the projected plane are 
computed as: 



^proj 



Jproj 



-X sin + y cos (f 

T] - 206265 Xpioj jr 
-X cos (/J sin - Y sin (/? sin + Z cos 9 

^ = 206265 yproj/r 



(6) 



Here 77 and ^ define the tangential coordinates on the sky in 
arcsec which are used to compute synthetic images comparable 
directly to observations as demonstrated in Section 6.2. 

The radial velocity and the position on the line of sight are 
given by: 



Vr - vx cos if cos + vy sin if cos 6 + vz sin 6 
Zi - X cos If cos + y sin cos + Z sin 



(7) 
(8) 



Then, 



- A grid corresponding to the desired map size and sampling 
is created. 

- In a single loop over all particles, we compute, which bin 
each particle will contribute to, and take into account this 
contribution. 

- In a second loop over the bins (pixels) we compute (if nec- 
essary) the actual values of a physical parameter. 

The algorithm is sufficiently fast to perform the map compu- 
tation in real time, e.g. the generation of a 400 x 400 pixels map 
takes a fraction of a second. In Fig. |6] we show an example of 
projected maps for one of the interactions. 



5.2. Spectrophotometric Properties 

We de veloped a technique using PEGASE.HR (iLe Borgne et aP 
120041) and PEGASE.2 dFioc & Rocca-Volmerangd Il997h to 



model spectrophotometric properties of interacting galaxies 
from the results of GalMer simulations. By taking into ac- 
count kinematics, star formation (SFH) and chemical enrich- 
ment history (CEH) we model spectra and broadband photomet- 
ric colours. We will present all the details regarding the spec- 
trophotometric modelling in a separate paper (Chilingarian et al. 
in prep.), but here we will briefly introduce the algorithms and 
results which can be obtained. 

There were several successful attempts 

dChakrabarti & WhitnevI 120091: iJonsson et all |20IO|) of mod- 
elling the spectrophotometric properties of the results of 
A^-body simulations. However, in all the known cases the stellar 
population information is taken into account in an approximate 
way, by characterising A^-body particles by their mean ages 
and metallicities. At the same time, the behaviour of spectral 
features is a strongly non-linear function of age and metallicity 
and it also differs significantly along the wavelength domain, in 
particular, due to the dust attenuation. Hence, SEDs computed 
using mass-weighted average ages and metallicities may not 
reflect real spectral energy distribution in galaxies. With the 
GalMer simulations, thanks to the usage of hybrid particles we 
are (I) able to trace SFH and CEH in detail through the entire 
duration of the simulation for every particle of this type; and 
(2) make use of this information. Therefore, we are able to 
make qualitatively better modelling of the spectrophotometric 
properties of interacting galaxies (and results of other A^-body 
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Fig. 6. Projected maps of 8 quantities traced by the GalMer sim- 
ulations for a merger of gSa and gSb, prograde orbit of type 5, 
650 Myr after the start of the simulations. Left-to-right, top-to- 
bottom: radial velocity, velocity dispersion {star and hybrid par- 
ticles), projected density of gas, projected density of old stars, 
metallicity of gas, metallicity of old stars, metallicity of new 
stars, total projected density {star + hybrid particles). 



TreeSPH simulations as well) and make direct comparison 
with observations. The SFH and CEH are computed using a 
whole sequence of snapshots in a given GalMer experiment as 
differences of gas masses in each hybrid particle in each of the 
10 metalUcity bins from [Fe/H]= -2.5 to +1.0 dex. 

The most time-consuming part of the modelling is running 
the PEGASE.HR code. Due to the complexity of the evolution- 
ary synthesis even with the present state of the art computer 
hardware it takes several seconds per spectrum (i.e. per star 
particle), resulting in several hours per snapshot including up- 
to 160000 hybrid and star particles. Computation of spectra for 
individual particles is absolutely crucial to properly model ef- 
fects of dust extinction and intrinsic broadening of spectral lines 
caused by motions of particles along a line of sight (i.e. inter- 
nal kinematics of galaxies). However, we can avoid the actual 
execution of the evolutionary synthesis code. 

The algorithms of evolutionary synthesis such as 
PEGASE.HR include: 

1. computation of isochrones for different ages and metallici- 
ties based on a given set of stellar evolutionary tracks and 
stellar IMF; 

2. picking up stellar spectra from a stellar library (either em- 
pirical or theoretical) for the atmosphere parameters corre- 
sponding to a given point on the isochrone; 

3. co-adding contributions of different types of stars according 
to the weights on the isochrone, making up simple stellar 



population (SSP) spectrum including stars of a single age 
and metallicity; 

4. co-adding different SSPs to reproduce complex SFH and 
CEH. 

For a case when the IMF is fixed and SFH and CEH are 
traced on a pre-defined grid of ages and metallicities, it becomes 
possible to execute only the 4th step from the list above. This 
means that if we pre-comp ute only once a grid o f SSPs corre- 
sponding to a given IMF ( iMiller & Scald (1 19791) in our case) 
and the grid of ages and metalUcities, we can reuse it for all the 
particles, improving the efficiency of the computations by sev- 
eral orders of magnitude. Once for every GalMer experiment 
we also pre-compute the SFH and CEH and store them as a two- 
dimensional histogram for every particle. 

For every spatial bin (see definition in the previous section) 
we first sort the particles along a line of sight to be able to ac- 
count for dust extinction. Then, either a high-resolution spec- 
trum or a broad-band spectral energy distribution (depending 
on the mode of the computation) is computed for every parti- 
cle starting from the most distant one from the observer. 

This is done by co-adding the pre-computed SSPs from the 
grid mentioned above with the weights corresponding to a mass 
contribution of stars of each age and metallicity contained in the 
SFH and CEH. 

The dust extinction is then applied to the total spectrum or 
multi-colour spectral energy distribution (SED) as it was com- 
puted at the previous step (i.e. excluding the current particle). 
For a known gas column density and a solar metallicity, we as- 
sume a standard dust-to-gas mass ratio to compute extinction as 
Av = A^h/1.871 X 10^' (Bohlinet al. 1978), where A^h is a num- 
ber of hydrogen atoms per cm', assuming Ry = Av/E{B - V) = 
3.1. We compute Nh along the line of sight and scale this for- 
mula linearly with metallicity. Then, we ac count for the wave - 
length dependence of extinction according to'Fitzpat rickl ( Il999l) . 
and apply it to the total stellar population SED generated at this 
step. 

After having applied the dust extinction, the total spectrum 
of a current particle is blue- or redshifted according to its radial 
velocity, which is done as a simple shift operation with linear 
interpolation, since our SSP grid is rebinned with a logarith- 
mic step in the wavelength corresponding to a fixed pixel size 
in km s"'. 

At the end, the total spectrum (or SED) of a current particle 
is co-added to the result. 

This algorithm makes it possible to compute in a few seconds 
for a given snapshot a total intermediate-resolution {R - 3000, 
3900 < A < 6800 A) spectrum, based on the ELODIE.3.1 
(iPrugniel et al.ll2007h empirical stellar library. The total number 
of co-added SSP in this extreme case reaches 10^. The broad- 
band F UV to NIR SED based on the low-resolution theoretical 
BaSeL (iLeieune et al.l 0*997) stellar library is computed much 
faster, because the number of SED points in it is only about a 
dozen compared to several thousands, allowing us to compute a 
400 X 400 pixels multi-colour datacube in a few seconds. 

Due to the limited spatial resolution of GalMer simulations, 
we are unable to reach the spatial scales sufficient to take into 
account properly the gas dumpiness, therefore we adopt a scale 
of 250 pc for the computation of extinction, which roughly cor- 
responds to the resolution of the simulations, and which results 
in the values of total extinction in simulated GalMer isolated 
galaxies well resembling observations. Our limited resolution 
also causes overestimation of extinction effects in central re- 
gions with high gas densities, where large amount of gas falls 
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Fig. 7. Examples of synthetic spectra generated by the spec- 
trophotometric modelling algorithm applied to the results of 
GalMer simulations. Two curves correspond to two different po- 
sitions along the projected major axis of a merger remnant: cen- 
tre approaching (blue), and receding (red). The spectra are nor- 
maUsed to unity at A = 5000 A for clarity. 



on-to during the interaction. Presently, we do not include any 
dust emission in mid- and far-IR spectral bands. This will be 
done in a future release of the GalMer database and will be pre- 
sented in a separate paper (Melchior et al. in prep.) 

In Fig. |7] we demonstrate an example of generated spectra 
for a merger remnant, where differences caused by the internal 
kinematics are clearly visible. Red and blue spectra taken 1 .5 kpc 
from the galaxy centre along the projected major axis, demon- 
strate the Doppler shift of spectral lines. Absorption lines are 
thus broadened naturally and reflect the intrinsic stellar velocity 
dispersion of the galaxy. 

At present, we do not attempt to model nebular emission 
lines, because our simulations do not trace at sufficient level 
of detail the physical conditions in the ISM. Only qualitative 
modelling can in principle be performed, assuming fixed ISM 
temperature and density depending only on the gas metallicity 
which is available for hybrid particles. This modelling will be 
addressed in detail in a separate paper (Melchior et al. in prep.) 



6. Astrophysical applications 

6. 1 . Galaxy properties from modelling 

At present, we explored the GalMer database to study a vari- 
ety of physical processes related to galaxy interactions, such as 
induced star formation enhancements, evolution of metallicity 
gradients, angular momentum redistribution, and its impact on 
the fin al kinematical propertie s of t he merger remnant. 

In iDi Matteo et"al1 (|2007|) and iDi Matteo et"al1 (l2008ah we 
investigated the enhancement of star formati on efficiency in 
galaxy interactions and mergers. We showed in lDi Matteo et al.l 



Lxy 

( l2007h that, in the final merging phase, retrograde encounters 
have greater star formation efficiency than direct encounters, that 
the amount of gas available in the galaxy is not the main param- 
eter governing the star formation efficiency in the burst phase, 
and that there is a negative correlation between the amplitude of 
the star forming burst at the merging phase and the tidal forces 
exerted at peric e ntral p assage. The general result presented in 
iDi Matteo et al.l (12008 ah shows that, at low redshift, galaxy in- 
teractions and mergers, in general, trigger only moderate star 
formation enhancements. Strong starbursts where the star for- 
mation rate is increased by a factor greater than 5 are rare and 
found only in about 15% of major galaxy interactions and merg- 
ers. Merger-driven starbursts are also rather short-lived, with a 
typical duration of activity of a few 10** yr. These conclusions 



are found to be robust, independent of the numerical techniques 
and star formation models. At higher redshifts, where galaxies 
are gas-rich, gas inflow induced starbursts are neither stronger 
nor longer than their local counterparts. These res ults are in good 
agreement with a number of observational works ( Bergvall et alj 
Uml lLietal.1120081: Ijogee^etaD 120091: Iknaoen & Jamesll2009h . 
demonstrating that interactions and mergers do not always trig- 
ger strong bursts o f star formation. 

More recently, dPi Matteo et a lJ2009bh . we investigated how 
the metallicity gradients in dry merger remnants depend on the 
structure and metallicity gradients of the galaxies involved in the 
merger. Our aim was to understand if dry mergers could lead to 
metallicity gradients as observed in elliptical galaxies in the lo- 
cal Universe, and if they always lead to a flattening of the initial 
(i.e., pre-merger) gradient. The analysis of the whole set of dry 
merger simulations in the GalMer database allowed us to show 
that the ratio of the remnant and the initial galaxy slopes spans 
a wide range of values, up to > 1 (with values greater than one 
resulting only when companions have gradients twice that of the 
progenitor). For a merger between two ellipticals having iden- 
tical initial metallicity slopes (i.e., equal companion and galaxy 
slopes), the metallicity profile of the remnant flattens, with a final 
gradient about 0.6 times the initial one. This flattening depends 
neither on the characteristics of the orbit of the progenitors or on 
their initial concentration. If the companion slope is sufficiently 
steep, ellipticals can maintain their original pre-merger metal- 
licity gradient. These results, compared to the observed variety 
of metallicity gradi ents in dwarf elliptical and lenticular galaxies 
(IChilin garianl2009"l) . may suggest the mergers to be an important 
channel of dE/dSO evolution. 

Given the diversity in outcomes of the mergers, we con- 
cluded that dry mergers do not violate any observational con- 
straints on the system atic characteristic s of metallicity gradients 
in local elHpticals (O gando et"al1l2005l) . 

The redistribution of the orbital angular momentum into 
internal rotation, and its impact on the kinematical proper- 
ties of the merger rernnant, b ecame a subject of two studies 
dDLMatteo et al."2008b','2009d). 



In Di Matteo et al. (2008bi), we presented a new scenario to 
form counter-rotating central components in early-type galaxies, 
by dissipative and dissipationless "mixed" mergers, consisting of 
elliptical-spiral systems in retrograde orbits. We demonstrated 
that the counter-rotation can appear both in dissipative and dis- 
sipationless retrograde mergers, and it is mostly associated to the 
presence of a disc component, which preserves part of its initial 
spin. In turn, the external regions of the two interacting galaxies 
acquire part of the orbital angular momentum, due to the action 
of tidal forces exerted on each galaxy by the companion. In the 
case of dissipative mergers, the central decoupled core could be 
composed of two distinct populations: the old stellar population, 
which has preserved part of its initial spin, and a new stellar 
population, born in situ from the kinematically decoupled gas 
component. 

Even the merger of two initially non-rotating, pressure sup- 
ported progenitor galaxies can lea d to remnant gal a xies ha ving 
peculiar kinematical properties. In iDi Matteo et al.l (l2009ah . we 
have indeed shown that it is possible to generate elliptical-like 
galaxies, with v/cr > 1 outside one effective radius, as a result of 
the conversion of orbital into internal angular momentum. This 
conversion occurs "outside-in": the external regions acquiring 
part of the angular momentum first, and it affects both the bary- 
onic and the dark matter components of the remnant galaxy (i.e. 
both acquire part of the angular momentum, the relative frac- 
tions depend on the initial concentration of the merging). If the 
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Fig. 9. Colour-magnitude (left) and colour-colour relations 
demonstrating contour density plots for 80000 observed low- 
redshift (z < 0.3) SDSS galaxies, GalMer galaxy models cor- 
responding to isolated galaxies of different masses and morpho- 
logical types, and a sequence of a GalMer simulation of a giant- 
giant interaction followed from the merger epoch for 1.5 Gyr 
every 50 Myr (red diamonds). Colours (see right panel legend) 
code morphological type, while symbol shapes (asterisk, trian- 
gle and square, see left panel legend) code luminosity-class. 




Fig. 8. RGB false-colour image of a simulated galaxy merger 
displayed over the SDSS image of the Abell 85 galaxy cluster. 



initial baryonic component is sufficiently dense and/or the en- 
counter takes place on a orbit with high angular momentum, 
the remnant galaxy exhibits hybrid properties, i.e. an elliptical- 
like morphology, but rotational support in the outer stellar halo 
(v/cr > 1). Systems with these properties have been recently 
observed through a combination of stellar absorption lines and 
planetary nebulae fo r kinematic studies of early-type galaxies 
(ICoccatoet all 120091) . Our results are in qualitative agreement 
with such observations and demonstrate that even mergers com- 
posed of non-rotating, pressure-supported progenitor galaxies 
can produce early-type galaxies with significant rotation at large 
radii. 

6.2. Synthetic observations: virtual telescope 

Some of the value-added data analysis services we provide for 
the GalMer database can be used as a "virtual telescope" to cre- 
ate simulated images and spectra of interacting galaxies. Thanks 
to the high-quality stellar population modelling provided by the 
PEGASE.2/PEGASE.HR code, we are able to compare the sim- 
ulated data to real observations. 

In Fig.[8j we present an RGB false-colour composite image 
of an interacting galaxy pair created from the results of spec- 
trophotometric modelling described above applied to GalMer 
simulations. The simulation result is superimposed over the 
ba ckground showing the Sloan Digital Sky Survey Data Release 
7 (lAbazaiian et al.ll2009l) image of the galaxy cluster Abell 85. 
We used SDSS g, r, and / bands in both observational and sim- 
ulated data, and the RGB visual isation code implementing the 
algorithm by iLupton et al.l (l2004l) to generate a false-colour im- 
age. Since the results of our spectrophotometric modelling are 
expressed in physical units, they can be transformed into ob- 
servables (mag arcsec"^) and projected on-to the sky (see eq|6]l, 
providing a direct way of comparing them to observational data. 
By co-adding model broadband fluxes and those coming from 



SDSS direct imaging, we can create realistically looking false- 
colour images such as the one shown in Fig.|8] 

Simulated images can then be analysed using classical ob- 
servational techniques, e.g. surface photom etry. A successful 
exampl e of such analysis is demonstrated in IChilingarian et al.l 
(I2009bh . where we found a match between a rather com- 
plex three-component observed density profile of the lenticular 
galaxy NGC 6340 and a merger remnant from the GalMer simu- 
lations, supporting the major merger to be an important event in 
the evolution of this particular object. 

We are also able to simulate the whole data cubes corre- 
sponding to the data produced by modern integral-field unit 
(IFU) spectrographs, such as SAURON at 4.2 m William 
Herschel Telescope, VIMOS at ESO Very Large Telescope, and 
GMOS at Gemini. Our spectral resolution and coverage allows 
us to model high-resolution blue setting of VIMOS (HR-Blue, 
R = 2500) and B600 grating of GMOS {R = 2200) usually cho- 
sen to study stellar kinematics of nearby galaxies. 

The spatial resolution of our simula tions (0.20-0.28 kpc) 
is comparable to that of the SAURON jde Zeeuw et al.ll2002l) 
survey targeting nearby galaxies at distances between 10 and 
40 Mpc, thus, having spatial resolution of 0.05 to 0.2 kpc given 
the average atmosphere seeing quality of 1 arcsec. At the same 
time, we can degrade the spatial resolution of our simulation 
in order to match observations made with VIMOS, e.g. post- 
starburst galaxies presented in Chilingarian et al. (2009a). 

As the stellar population models improve, we will be able to 
upgrade the spectrophotometric modelling engine in the GalMer 
database in order to simulate data coming from next genera- 
tion IFU facilities such as the second-generation VLT instrument 
MUSE. 

6.3. Colour-magnitude relations 

The spectrophotometric modelling facilities allow us to compute 
total magnitudes and colour of interacting galaxies and to di- 
rectly compare them to observational data. 

The colour-magnitude and colour-colour plots comparing 
results of GalMer simulations to a large observational dataset are 
pre sented in Fig. [9] A we ll-known galaxy colour bimodality (see 
e.g. IStrateva et aLfcooll) is clearly seen in the contour plot pre- 
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senting the distribution of some 80000 nearby galaxies (z < 0.3) 
from the SDSS survey. The photometric data were corrected for 
the Galactic extinction according to Schlegel et al. (1998!) and 
converted into rest-frame, i.e. A:-corrected using the analytical 
approximations presented in Chilingarian et al. (submitted). The 
"red sequence" and the "blue cloud" separ ated by the "green val- 
ley", a locus of post-starburst galaxies (e.g. lGoto et alj2003l) can 
be identified. 

The total magnitudes of GalMer mergers in SDSS bands 
are computed using the spectrophotometric modelling described 
above. We show the positions of GalMer galaxy models (isolated 
galaxies) of different masses and morphological types. Elliptical 
and lenticular galaxies (red and yellow symbols) are sitting on 
the red sequence, while spiral galaxies are bluer, directed down 
toward the blue cloud. GalMer Sa and Sb spiral galaxies (green 
and blue symbols) at an evolutionary stage shown in Fig. preside 
in the green valley. 

We also follow the evolution of an interacting pair from 
the time when galaxies merged (red diamonds connected with 
lines) till the end of the simulation showing intermediate snap- 
shots every 50 Myr. Right after the gas-rich merger, remnants 
sit in the "blue cloud", despite very strong internal extinction of 
regions of massive star formation. Then, as the star formation 
decreases, luminosity-weighted ages of stellar populations in- 
crease causing global colours to become redder and, thus, mov- 
ing merger remnants in the top-right direction across the "blue 
cloud" through the "green valley" toward the "red sequence". 
Finally, 400-600 Myr after the merger time, all remnants settle 
on-to the "red sequence", although some residual star formation 
is still observed. 

One can see that merger remnants cross the "green valley" 
very fast giving possible explanation why there is at present a 
deficit of galaxies in this region of the colour-magnitude dia- 
gram, given that all mergers in this mass range have not yet been 
included in the database. 

7. Summary 

We present the GalMer database providing online access to the 
results of TreeSPH simulations. The structure of the database 
conforms to recent International Virtual Observatory standards. 
We describe the interactive data access web-interface, advanced 
mechanisms for data visualisation and manipulation using con- 
nection to the Virtual Observatory tools dedicated for dealing 
with tabular, imaging and spectral data. 

After having stored the snapshots in FITS binary table for- 
mat, giving access to all coordinates and data on individual sim- 
ulation particles, we provide a set of value-added tools using the 
results of TreeSPH simulations. These include; (1) generator of 
projected maps of various physical quantities traced by the simu- 
lations; (2) engine to perform the modelling of spectrophotomet- 
ric properties of interacting and merging galaxies based on the 
PEGASE.2/PEGASE.HR evolutionary synthesis code. The lat- 
ter tool can be used as a virtual telescope and synthetic images, 
spectra, and datacubes generated using it are directly compara- 
ble to observational data. 

We provide examples of several use-cases for the GalMer 
database using our value-added tools. 

The database will be updated by including new simulations 
as they have been completed. 
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Appendix A: Galaxy models: density profiles 



In Sect l2.1l we presented the galaxy models adopted for the sim- 
ulations. In particular, we saw that the dark halo and the o ptional 
bulge are modelled as a Plummer sphere (Binney & Tr emainej 
Il987l pag.42), with characteristic masses Mb and Mh and char- 
acteristic radii and r//. Their densities are given, respectively, 
by the following analytical formula: 



pH(r) 



and 



.2 \-5/2 



.2 \-5/2 



(A.l) 



(A.2) 



The ste llar and gaseous discs follo w a Miyamoto-Nagai density 
profile (iBinnev & Tremaindl 19871 pag.44): 



(A.3) 



An 



QgR- + {ag + 3 ^z^ + hg^) {ug + yjz^ + hg^^ 



(A.4) 



with masses M, and Mg and vertical and radial scale lengths 
given, respectively, by /i, and «», and hg and Og. 



Appendix B: Galaxy models evolved in isolation 

Since the aim of this work is to study the effects that interactions 
and mergers have in driving galaxy evolution, it is important to 
study, for comparison, also the evolution of the galaxies in iso- 
lation, in order to distinguish properly secular processes from 
those related to the interaction. In Figs. IB. II and IB. 21 we show 
the evolution of the gaseous and stellar components of the gi- 
ant gSa, gSb and gSd galaxies, respectively, evolved in isolation 
for 3 Gyr As it can be seen, all the models evolve rapidly in the 
first 0.5-1 Gyr; a stellar bar and spiral arms form, and at the same 
time gas compresses into density waves and clumps and partially 
falls into the central regions, where a star formation enhance- 
ment takes place (fo r the star formation ra te of these galaxies, 
we refer the reader to lDi Matteo et al.ll2007h . 
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Fig.B.l. Gas maps of the gSa, gSb, and gSd models with Qgas - 
0.3 evolved in isolation. Snapshots are equally spaced in time 
from t=0 to t=3 Gyr. Each map is 60 kpc x 60 kpc in size. 




Fig. B.2. Star maps of the gSa, gSb, and gSd models evolved in 
isolation. Snapshots are equally spaced in time from t=0 to t=3 
Gyr. Each map is 60 kpc x 60 kpc in size. 



